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cyclically  shifted,  then  truncated  and  transmitted.  Although  shortened,  the  trucated  code  still  represents  k  bits  of  information,  thus  leading  to  an 
increased  data  rate.  We  evaluate  radiometer  detection  of  TCCSK  and  it  is  shown  that  the  probability  of  detection  is  increased  compared  with  the 
detection  of  CCSK. 
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A  low  probability  of  intercept  (LPI),  or  low  probability 
of  detection  (LPD)  communication  technique  known  as  cyclic 
code  shift  keying  (CCSK)  is  described.  We  discuss  the  basic 
concepts  of  CCSK  and  describe  a  system  based  on  the  use  of 
random  or  pseudorandom  codes  for  biphase  modulation.  We  use 
simulation  to  show  that  the  bit  error  rate  (BER)  for  CCSK  can  be 
closely  estimated  by  using  existing  equations  that  apply  to  M-ary 
orthogonal  signaling  (MOS).  Also,  we  show  that  significantly 
fewer  computations  are  required  for  CCSK  than  for  MOS 
when  the  number  of  bits  per  symbol  is  the  same.  We  show  that 
using  biphase  modulation  results  in  waveforms  that  have  a  large 
time-bandwidth  product  and  very  low  input  signal-to-noise  ratio 
(SNR)  and  thus  inherently  have  an  LPI  by  a  radiometer.  We 
evaluate  detection  by  a  radiometer  and  show  that  LPI  can  be 
achieved  by  using  codes  of  lengths  greater  than  about  212  (i.e., 
by  transmitting  more  than  about  12  bits  per  symbol).  Results 
illustrate  the  effect  that  the  CCSK  symbol  length  and  error 
probability,  and  the  radiometer  integration  time  and  probability 
of  false  alarm  (PFA),  have  on  detection  by  a  radiometer.  We 
describe  a  variation  of  CCSK  called  truncated  CCSK  (TCCSK). 

In  this  system,  the  code  of  length  2*  is  cyclically  shifted,  then 
truncated  and  transmitted.  Although  shortened,  the  truncated 
code  still  represents  k  bits  of  information,  thus  leading  to  an 
increased  data  rate.  We  evaluate  radiometer  detection  of  TCCSK 
and  it  is  shown  that  the  probability  of  detection  is  increased 
compared  with  the  detection  of  CCSK. 
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Cyclic  code  shift  keying  (CCSK)  is  a  form  of 
M-ary  signaling  over  a  communication  channel  [1], 

In  its  simplest  form,  a  “base  function”  f(t )  is  chosen, 
and  a  cyclically  (circularly)  shifted  version  of  / ( t )  is 
used  to  modulate  a  carrier.  The  function  /(f)  has  the 
property  that  its  cyclic  autocorrelation  has  a  distinct 
peak  and  “low”  sidelobes.  Assuming  synchronization, 
the  receiver  cyclically  correlates  the  received  signal 
plus  noise  with  /(f)  and  estimates  the  position  of 
the  correlation  peak.  If  the  number  of  resolvable 
positions  is  M,  the  number  of  bits  per  “symbol”  is 
B  -  log 2M.  The  base  function  we  consider  here  is  a 
binary  sequence  b  =  (b0,bx,...,bM_x)T  with  bm  =  ±1, 
resulting  in  biphase  modulation,  of  the  carrier.  We 
also  describe  a  technique  that  increases  the  bit  rate  by 
using  truncations  of  shifted  versions  of  b  to  biphase 
modulate  the  carrier.  This  technique  is  referred  to  as 
truncated  cyclic  code  shift  keying  (TCCSK). 

Three  methods  for  generating  b  are  discussed, 
including  a  maximal-length  sequence  (MLS)  [2],  a 
modified  maximal-length  sequence  (MMLS)  [3]  and  a 
randomly  chosen  sequence.  Maximal-length  sequence 
generators  produce  sequences  with  elements  + 1  and 
0,  which  we  convert  to  ±1  by  replacing  the  zeros  with 
—  1.  However,  we  use  the  term  MT  S  to  describe  both 
types  of  sequences  when  the  type  is  clear  from  the 
context. 

Conventional  M-ary  orthogonal  signaling  (MOS) 
[4,  p.  167  ff.]  uses  one  of  M  =  2k  orthogonal 
functions  to  modulate  a  carrier.  The  receiver  correlates 
the  received  signal  plus  noise  with  each  of  the 
orthogonal  functions  and  determines  the  one  with 
the  highest  correlation.  We  show  by  Monte  Carlo 
simulation  that  CCSK  performance  in  Gaussian  noise 
is  essentially  the  same  as  MOS,  when  the  symbol 
error  probability  Ps  is  larger  than  about  10~4. 

The  primary  reason  for  choosing  CCSK  for 
M-ary  signaling  instead  of  conventional  MOS  is 
the  simplicity  of  the  signal  processing.  We  show 
that  CCSK  only  requires  the  computation  of  the 
Fourier  transform  of  the  received  signal  plus  noise 
followed  by  an  inverse  transform  of  the  product 
of  this  transform  and  the  complex  conjugate  of 
the  transform  of  the  base  function.  For  the  codes 
considered  here,  these  operations  are  performed  by 
using  the  discrete  Fourier  transform  (DFT)  or  the  fast 
Fourier  transform  (FFT).  This  is  contrasted  with  the 
processing  for  MOS,  which  requires  the  correlation 
of  each  of  the  M  orthogonal  functions  (e.g.,  Walsh 
functions)  with  the  received  signal  plus  noise. 

The  use  of  the  binary  sequence  b  described  above 
as  the  base  function  leads  to  a  system  that  provides 
a  low  probability  of  intercept  (LPI)  by  a  radiometer, 
because  of  the  large  processing  gain.  We  briefly 
discuss  the  characteristics  of  a  radiometer  and  derive 
the  basic  equations  used  to  evaluate  its  performance. 
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These  equations  are  then  applied  to  the  CCSK  and 
TCCSK  systems  to  evaluate  their  vulnerability  to 
detection. 


II.  CYCLIC  CODE  SHIFT  KEYING 

CCSK  uses  discrete  cyclic  shifts  of  a  base  function 
/(f)  to  modulate  a  carrier,  with  each  shift  representing 
B  bits  of  information.  The  receiver  determines  the 
shift  by  computing  the  cyclic  correlation  of  the  base 
function  /(f)  with  the  received  signal  plus  noise.  The 
base  function  we  consider  here  is  a  binary  sequence 

b  =  (b0,b1 . bu_ i)r  with  bm  =  ±1,  resulting  in 

biphase  modulation  of  the  carrier.  A  typical  sequence 
b  is  shown  in  Fig.  1 .  The  elements  of  b  indicated  in 
Fig.  1  are  called  “chips.”  Each  chip  is  of  duration  Tc 
and  the  “symbol”  duration  Ts  is  MTC.  We  define  the 
chip  bandwidth  as  Wc  -  l/Tc,  so  that  TCWC  —  1. 

We  consider  three  different  methods  for  generating 
b.  The  first  uses  an  MLS  of  length  M  =  2k  -  1,  which 
has  the  property  that  it3  cyclic  autocorrelation  has  a 
peak  of  M  and  sidelobes  of  - 1 .  Unfortunately,  in 
this  case  the  number  of  bits  B  is  less  than  k.  Also, 

M  is  not  a  power  of  two,  thus  processing  by  using 
the  FFT  algorithm  may  be  precluded.  To  alleviate 
these  problems,  an  MMLS  (also  called  an  “extended 
m-sequence”)  [3]  is  used.  In  this  case  an  MLS  is 
generated  and  a  - 1  or  + 1  is  inserted  to  extend  the 
length  to  M  =  2*.  A  result  of  this  modification  is  an 
increase  in  the  level  of  the  autocorrelation  sidelobes, 
compared  with  the  true  MLS.  We  show  later  that  this 
increase  has  little  effect  on  error  probabilities  because 
the  level  of  the  input  signals  of  interest  is  so  low 
that  the  error  process  is  controlled  by  the  noise.  This 
fact  leads  us  to  consider  a  third  option  in  which  b  is 
obtained  by  generating  a  random  binary  sequence  of 
±1  s. 


A.  Noncoherent  CCSK  Using  Binary  Sequences 

The  cyclic  shifts  of  the  code  b  are  designated 
b0,b1,...,bAjf_I,  where  b0  =  b,  and  b„  is  the  nth  shift. 
We  consider  the  case  where  it  is  reasonable  to  assume 
time  synchronization  so  that  the  receiver  knows  the 
time  of  arrival  of  each  “symbol”  to  a  small  fraction 
of  one  chip  interval.  We  also  assume  that  the  carrier 
frequency  is  known.  However,  the  relative  phase  of 
the  received  signal  is  not  known,  and  noncoherent 
processing  is  required.  We  assume  the  received  signal 
is  a  sum  of  signals  of  the  form 

r(t)  =  A  cos(uj0(  +  <pmn  +  0)  +  g(t)  (I) 


where  A  is  the  amplitude  of  the  received  signal, 
ipmn  =  0  or  7r  is  determined  by  the  mth  element  of 
the  transmitted  code  bn,  9  is  the  unknown  phase,  and 
g(t)  is  a  Gaussian  noise  process  with  mean  zero  and 
(one-sided)  power  spectral  density  N0/2.  A  quadrature 
(or  baseband)  detector  is  used  and  its  sampled  output 
is  a  column  vector 

r  =  Ab„exp|j0]  +g.  (2) 

In  (2),  bn  is  a  column  vector  of  the  shifted  code 
elements,  g  is  a  column  vector  of  independent  and 
identically  (HD)  circular  Gaussian  noise  with  variance 
a2,  and  A  is  the  signal  amplitude.  The  receiver 
computes 

-  |bmr|  =  |Ab£b„  +  b^g|  (3) 

for  m  -  0, 1, . . . ,M  -  1 .  If  max{Ym}  =  Sp,  the  bits 
corresponding  to  the  pth  symbol  are  output.  A  symbol 
error  occurs  if  p  /  n. 

B.  Comparisons  with  M- ary  Orthogonal  Signaling 

In  conventional  MOS,  one  of  a  set  of  M 
orthogonal,  equal-energy  signals  is  used  to  represent 
B  =  log2M  bits  of  information.  One  form  of  the 
receiver  correlates  the  received  signal  plus  noise  with 
each  of  the  M  reference  signals  and  determines  which 
has  the  maximum  correlation.  Although  a  single  cyclic 
correlation  is  performed  in  implementing  CCSK, 
the  process  can  be  viewed  as  the  performance  of  M 
separate  correlations,  as  implied  in  (3).  If  the  cyclic 
shifts  of  the  code  were  all  uncorrelated,  then  CCSK 
would  be  a  special  case  of  MOS.  However,  for  codes 
that  provide  LPI,  the  cyclic  shifts  are  correlated;  that 
is,  the  cyclic  autocorrelation  function  has  non-zero 
“sidelobes.”  We  show  that  the  bit  error  rate  (BER)  for 
CCSK  can  be  approximated  by  using  the  equations  for 
MOS,  even  though  the  cyclic  shifts  are  correlated. 

C.  Performance  Estimates 

The  performance  of  noncoherent  MOS  in  terms  of 
BER  is  well  documented  in  the  literature  [5,  p.  489]. 
For  M  orthogonal  signals,  the  probability  of  a  symbol 
(or  word)  error  is  given  by 

poo 

Ps  =  1  —  /  exp  [-(x+q)]I^(-\/4qx)H(x)dx  (4a) 

Jo 

where 

H(x)  =  [1  -  (4b) 

q  =  Es/Nq  is  the  ratio  of  symbol  energy  to  noise 
power  spectral  density,  and  I0(x)  is  the  modified 
Bessel  function  of  the  first  kind  and  order  zero.  By 
using  the  binomial  expansion  of  H ( x )  and  integrating 
(4a)  term-by-term,  a  finite  series  for  Ps  is  obtained 
[5,  p.  489].  However,  the  series  is  alternating  in  sign 
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Fig.  2.  Comparison  of  CCSK  performance  with  MOS,  M  =  1024. 


Fig.  3.  Comparison  of  CCSK  performance  with  MOS,  M  =  4096. 


and  presents  numerical  difficulties  for  large  M  or  q. 
Fortunately,  (4a)  lends  itself  to  numerical  integration, 
and  that  is  the  technique  we  use  for  its  evaluation. 
With  M  =  2k,  the  probability  PB  of  a  bit  error  is 
obtained  directly  from  Ps  as  [5,  p.  198] 


2*-i 

pB-ps2rzr\- 


(5) 


Monte  Carlo  simulations  were  performed  to 
compare  CCSK  performance  with  MOS.  Some 
of  the  results  are  shown  in  Figs.  2-4.  From  these 
results  we  conclude  that  (4)  can  be  used  as  a  close 
approximation  to  the  symbol  error  probability 
for  CCSK  when  Ps  is  larger  than  about  10-4  and 
when  M  >  1024.  This  conclusion  is  based  on  the 
simulations  performed;  however,  there  appears  to  be 
no  reason  to  assume  that  the  approximation  is  not 


valid  for  smaller  Ps,  especially  for  large  values  of 
M.  To  corroborate  the  simulation  results,  we  show 
that  the  cross-correlation  between  cyclic  shifts  are 
small  and  thus  are  “approximately”  orthogonal. 

Table  I  shows  the  maximum  cross-correlation  for 
k  =  9, 10,. ..,16  for  MLS,  MMLS,  and  random 
sequences.  The  MMLS  results  were  obtained  by  using 
a  single  set  of  shift-register  taps  for  each  k  to  generate 
the  sequence.  Some  slight  variations  are  expected  . 
if  other  sets  of  taps  are  used.  Note  that  for  k  >  12, 
the  maximum  cross-correlation  is  less  than  0.03  for 
MMLS  and  decreases  with  increasing  k.  The  tabulated 
value  for  each  random  sequence  was  obtained  by 
averaging  the  results  from  100  different  random 
sequences.  A  comparison  shows  that  the  random 
sequence  has  larger  maximum  cross-correlation  than 
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Fig.  4.  Comparison  of  CCSK  performance  with  MOS,  M  =  16384. 


TABLE  I 

Maximum  Cross-Correlation  of  Cyclic  Shifts,  M  =  2k 


k 

9 

10 

11 

12 

■  13 

14 

15 

16 

MLS 

0.0002 

2-14  —  1 

2-15  _  i 

2-16  _j 

MMLS 

0.0625 

0.0586 

0.0391 

0.0254 

RANUUM 

0.0760 

0.0563 

0.0422 

6.0308 

0.0228 

the  MMLS.  As  expected,  the  MLS  has  the  least 
cross-correlation. 

The  simulation  results  presented  in  Figs.  2-4, 
along  with  the  data  in  Table  1 ,  give  credence  to 
concluding  that  the  error  probabilities  for  MOS  are 
good  estimates  of  those  for  CCSK.  This  is  especially 
true  for  larger  values  of  k  because,  as  indicated, 
the  maximum  cross-correlation  decreases  for  both 
MMLS  and  random  sequences.  Also,  this  estimate  is 
applicable  when  the  input  signal-to-noise  ratio  (SNR) 
is  low  and  the  additive  noise  obscures  the  non-zero 
cross-correlation  between  the  cyclic  shifts  of  b. 

D.  Implementation  Issues 

Equation  (3)  is  a  representation  of  the  processing 
required  by  the  receiver  and  is  useful  in  comparisons 
with  MOS.  However,  in  practice  we  compute  the 
cyclic  correlation  by  using  the  DFT.  The  DFT  of 
the  code  b  is  computed  and  its  complex  conjugate  is 
stored.  The  receiver  computes  the  DFT  of  the  sampled 
received  signal  r,  which  is  defined  in  (2),  and  obtains 

S  =  iIDFT(DFT*  (b)  x  DFT(r))  |  (6) 

where  IDFT  is  the  inverse  DFT.  (The  multiplication 
implied  in  (6)  is  term-by-term.)  The  elements  of  the 
vector  S  obtained  by  using  (6)  are  S0,S1,...,  Sm-i> 
which  are  the  same  as  are  obtained  by  using  (3).  From 


(6),  it  follows  that  the  processing  required  for  CCSK 
can  be  significantly  less  than  for  MOS,  which  requires 
the  computation  of  all  M  individual  cross-correlations. 

E.  TCCSK 

The  CCSK  system  previously  described  transmits 
all  M  elements  (chips)  of  cyclic  shifts  of  the  code  b  to 
represent  symbols.  Thus,  the  data  rate  is  1  /Ts  symbols 
per  second,  or  k/Ts  bits  per  second.  (See  Fig.  1.) 

We  may  increase  the  bit  rate  while  maintaining  the 
same  bandwidth  by  using  the  technique  we  refer 
to  as  TCCSK.  As  discussed  in  the  next  section,  the 
code  sequence  b  is  chosen  to  have  the  property  that 
its  cyclic  autocorrelation  has  a  distinct  peak  and  low 
sidelobes.  If  the  code  length  is  large  (e.g.,  M  >  1024), 
then  a  subsequence  br  of  b  exhibits  these  properties 
when  cyclically  correlated  with  b.  This  fact  leads  to 
the  use  of  TCCSK. 

Instead  of  transmitting  all  M  chips  of  the  shifted 
code  bn  to  represent  the  nth  symbol,  the  code  bn 
is  truncated  and  only  the  first  MT  chips  bnT  of  hn 
are  transmitted.  The  received  signal  plus  noise  is 
cyclically  correlated  with  b  by  first  appending  M  —MT 
zeroes  to  rr  =  AbnT  exp[/0]  +  gr.  (See  (2).)  As  a 
result,  the  number  of  bits  per  symbol  is  still  k  but  the 
bit  rate  has  been  increased  by  a  factor  of  M/MT.  If 
TCCSK  is  to  maintain  the  same  BER  as  CCSK,  then 
Es/N0  must  be  the  same  for  both.  This  means  that  the 
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amplitude  of  the  signal  must  be  increased  by  a  factor  by 
of  (M [MTyi2.  This  is  considered  later  when  detection 
by  a  radiometer  is  discussed. 


PMR  = 


M2A2+Mo2 
QA2  +  Mo2 


(10) 


III.  CODE  SELECTION  TRADEOFFS 


where  Q  is  the  average  sidelobe  power  of  the  code  b. 
Equation  (10)  can  be  written  as 


For  CCSK,  the  information  transmitted  is 
contained  in  the  location  of  the  maximum  of  the 
correlation  function  defined  by  (3)  and  (6).  Thus, 
a  measure  of  performance  that  can  be  used  in  code 
selection  is  the  peak-power  to  mean-sidelobe-power 
ratio  (PMR).  If  code  b'  has  a  significantly  higher 
PMR  than  code  b",  then  it  is  intuitive  that  b'  will 
provide  better  performance.  Another  consideration 
in  code  selection  is  their  noise-like  property.  MLSs  are 
sometimes  referred  to  as  “pseudonoise”  sequences, 
but  they  possess  certain  structure  not  found  in  a 
truly  random  sequence.  In  an  MLS  (with  elements 
1  and  0)  of  length  M  =  2k  -  1,  all  k-bit  binary 
numbers  except  zero  appear  as  k  successive  elements. 
For  example,  one  MLS  of  length  seven  is  b  = 
[101001  1].  Each  three  successive  bits  are 
binary  representations  of  [5  2  4  1  3  7  6],  with 
the  last  two  [7  6]  determined  cyclically.  Some  of 
this  same  structure  still  is  retained  if  an  MMLS 
is  used. 

A.  PMR 

To  simplify  the  analysis  in  defining  and  evaluating 
the  PMR  of  a  sequence  b,  we  assume  that  the  phase 
of  the  received  signal  is  known.  In  this  case,  the 
received  signal  plus  noise  is  r  =  Ahn  +  g,  where  g 
is  a  vector  of  HD  Gaussian  noise  with  variance  a2. 
Without  loss  of  generality,  we  may  assume  that  n  =  0; 
i.e.,  assume  that  bn  =  b.  The  correlation  of  r  with 
b  produces  M  terms.  The  first  term,  defined  as  the 
“peak,”  is 

s0  =  MA  +  h0  (7) 

where  h0  is  a  Gaussian  random  variable  with  mean 
zero  and  variance  Mu2.  The  remaining  M  —  1  terms 
(the  sidelobes)  are  of  the  form 


1  +MSNRin 
1  +  QSNR^/M 


(ID 


where  SNRM  =  A2/cr2  is  the  input  SNR. 

Fig.  5  shows  PMR  versus  SNRjn  for  an  MLS, 
MMLS,  and  a  random  binary  (RANDOM)  sequence. 
Note  that  the  MLS  curve  is  basically  a  unity-slope  line 
up  to  about  20  dB  input  SNR.  What  is  most  striking, 
however,  is  the  fact  that  the  three  curves  are  nearly 
coincident  for  negative  input  SNR,  which  is  the  usual 
case  for  use  in  CCSK.  Also,  it  is  obvious  from  the 
figure  that  a  unity-slope  line  provides  an  excellent 
estimate  of  PMR  for  the  region  of  interest  (negative 
input  SNR).  Further  results  in  the  next  figure  illustrate 
these  points. 

The  PMR  versus  SNRin  for  M  =  210,  213  and  216 
is  shown  in  Fig.  6  for  a  random  sequence  and  an 
MMLS.  Note  that  the  three  pairs  of  curves  follow 
the  unity-slope  line  up  to  SNR1n  of  about  -5  dB.  An 
exception  is  the  curve  for  210  when  SNR1n  is  less 
than  about  -20  dB,  the  region  where  the  position 
of  the  peak  correlation  is  likely  to  be  determined  by 
noise.  That  is,  the  peak  does  not  necessarily  occur  at 
n  =  0. 


B.  PMR  for  TCCSK 

When  TCCSK  is  used  it  is  necessary  to  modify 
(11)  to  compute  the  PMR.  Because  only  MT  chips  of 
the  code  are  transmitted,  the  peak  correlation  is  given 
by 

5q  j  —  MjA  +  k(yj-  (12) 

where  h0T  is  a  Gaussian  random  variable  with  mean 
zero  and  variance  MTcr2.  Similarly,  the  sidelobes  are  of 
the  form 

SmT  =  +  \ mT  (^) 


sm  +  hm 


(8) 


where  qm  is  the  mth  signal  sidelobe  and  hm  is  a 
Gaussian  random  variable  with  mean  zero  and 
variance  Mo2.  We  define  PMR  as 


PMR  = 


g[?g j 

AV{E[s2]} 


(9) 


where  AV  denotes  the  average  over  the  M  —  1 
sidelobes.  Note  that  we  define  the  peak  to  be  located 
at  the  position  of  the  shifted  code  (in  this  case, 
zero);  however,  for  low  SNR,  the  actual  maximum 
may  occur  at  some  other  position  (i.e.,  an  error 
occurs).  By  using  (7),  (8),  and  (9),  the  PMR  is  given 


where  qmT  is  the  mth  signal  sidelobe  and  hmT 
is  Gaussian  with  mean  zero  and  variance  MTo2. 
Therefore,  the  PMR  of  the  truncated  sequence  is 


PMRr  = 


1  +MrSNRIN 
1  v  QrS  NRin  /  A4-- 


(14) 


where  QT  is  the  average  sidelobe  power  of  the 
truncated  code. 

Fig.  7  shows  PMRr  as  a  function  of  the  input 
SNR  for  M  =  213  and  MT  =  213,  212,  and  2n.  (When 
Mt  =  213,  there  is  no  truncation.)  This  shows  that 
PMRr  decreases  by  about  3  dB  when  Mr  is  decreased 
by  a  factor  of  two.  However,  this  fact  is  misleading 
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Fig.  6.  PMR  versus  input  SNR  for  CCSK. 


when  error  performance  is  considered.  To  maintain 
the  same  error  performance,  EB/NQ  must  be  kept 
constant  when  MT  is  decreased;  that  is,  the  input  SNR 
must  increase.  Fig.  8  shows  PMRr  as  a  function  of 
Eb/N0  over  the  range  of  0  to  20  dB.  Note  that  for 
small  Eb/N0,  PMRr  is  essentially  the  same  for  both 
truncated  and  nontruncated  sequences.  In  fact,  as 
shown  in  the  next  section,  the  range  of  interest  for 
Eb/N0  is  about  0  to  5  dB,  and  over  this  range  PMR 
and  PMRr  differ  by  only  a  fraction  of  a  dB. 

C.  Random  versus  Pseudorandom  Codes 

By  observing  graphs  of  the  probability  of  a  bit 
error  PB  versus  Es/N0  (e.g.,  [5,  Fig.  10-6])  we  see  that 


for  10  <  k  <  20  and  0.00001  <  PB  <  0,1,  the  required 
Eb/N0  is  less  than  about  4.5  dB.  Also,  EB/N0  (in  dB) 
is  given  by 

Eb/N0  =  SNRin  +  101og10M  -  101og10*  (15) 

which  means  that  SNR]N  is  less  than  about  -15:6  dB. 
(For  TCCSK,  EB/NQ  is  given  by  (15)  with  M  replaced 
by  MT.)  The  data  in  Fig.  6  show  that  the  PMR  for 
MMLS  and  a  random  sequence  is  essentially  the 
same  over  the  range  of  SNRin  of  interest.  Therefore, 
random  sequences  can  be  used  as  the  code  in  CCSK 
with  little  effect  on  performance  compared  with 
MMLS. 

Some  caution  must  be  used  in  generating  a  random 
sequence  for  the  CCSK  code.  For  example,  although 
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the  sequence  is  called  “random,”  it  is  likely  to  be 
obtained  by  using  a  random  number  generator  that 
is  actually  “pseudorandom.”  Thus  the  quality  of  the 
generator  used  must  be  ensured  [6],  (The  random 
codes  used  in  Figs.  5  and  6  were  generated  by  using 
the  Matlab*  routine  “rand.”)  Additionally,  some 
random  codes  may  have  unacceptable  properties. 

For  example,  it  may  be  necessary  to  ensure  that  the 
numbers  of  1  s  and  —  1  s  are  approximately  equal  to 
avoid  having  a  dc  offset. 

IV.  VULNERABILITY  OF  CCSK  TO  RADIOMETER 
DETECTION 

A  radiometer  (or  energy  detector)  is  often  the  most 
effective  device  to  detect  spread-spectrum  signals  [7]. 


It  is  conceptually  a  simple  device,  and  requires  only 
a  few  assumptions  to  be  made  about  the  structure  of 
the  signals  being  detected.  Invulnerability  to  detection 
by  a  radiometer  is  required  if  a  communication  system 
is  to  be  considered  LPI.  Fig.  9  is  a  simplified  block 
diagram  of  a  radiometer. 

A.  Equations  for  Evaluating  Performance 

Equations  for  evaluating  radiometer  performance 
are  given  in  [7]  and  are  based  on  derivations  by 
Urkowitz  [8].  We  assume  that  the  noise  at  the 
input  to  the  radiometer  is  a  zero-mean,  stationary, 
Gaussian  random  process  that  has  a  flat,  bandlimited 
(one-sided)  power  spectral  density  iV0/2  over  the 
bandwidth  W  of  the  bandpass  filter.  For  convenience, 
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Fig.  9.  Radiometer  block  diagram. 


we  consider  a  normalized  detection  statistic 
V  =  2V /N0.  Urkowitz  shows  that  the  probability 
distribution  of  V'  is  closely  approximated  by  the 
chi-square  distribution  with  2 TW  deg  of  freedom 
in  the  noise-only  case  and  by  the  noncentral 
chi-square  distribution  with  2 TW  deg  of  freedom 
and  noncentrality  parameter  A  =  2ER/N0  =  2 S  in  the 
signal-plus-noise  case.  The  parameter  Ep  is  the  total 
signal  energy  integrated  by  the  radiometer  and  Ar0  is 
the  noise  power  spectral  density.  For  simplicity,  we 
assume  that  TW  is  an  integer. 

Detection  is  accomplished  by  comparing  the 
normalized  radiometer  output  V  with  a  threshold  K, 
and  a  signal  is  claimed  present  if  V'  >  K .  If  V'  >  K 
when  no  signal  is  present,  a  false  alarm  occurs.  By 
making  appropriate  changes  of  variables  in  the  central 
and  noncentral  chi-square  distributions,  and  defining 
Kq  =  K/2,  the  equations  for  the  probability  of  a  false 
alarm  PFA  and  the  probability  of  detection  PD  are 
given  by  [7,  pp.  56-57] 


PFA 


„TW- 


(TW-1)\ 


du 


(16) 


and 

(17) 

where  j  (x)  is  the  modified  Bessel  function  of  the 
first  kind  and  order  TW  —  1.  These  equations  were 
used  to  evaluate  the  detectability  of  the  CCSK  biphase 
modulated  waveforms  considered  here. 


B.  Radiometer  Detection  of  CCSK 

To  simplify  the  analysis,  we  assume  that  the 
radiometer  has  the  same  gains  and  losses  as  the 
intended  communications  receiver,  and  is  located  at 
the  same  range.  As  a  result  of  this  assumption  both 
have  the  same  SNRjN,  so  that  ER  =  Es  when  T  =  TS. 
We  also  assume  that  the  radiometer  integration  time 
varies  in  units  of  the  chip  time  Tc.  Then,  because 
TCWC  =  1,  the  radiometer  TW  product  is  just  the 
number  of  chips  integrated.  Also,  ER/N0  (in  dB)  can 
be  expressed  as 

Er/N0  =  SNRjn  +  101ogIO  TW.  (18) 

By  using  (15)  and  (18),  we  have  the  following  relation 
between  ER/N0  and  EB/N0: 

Er/N0  =  EB/N0  +  101og10TW  -  101og10M  +  101og10*. 

(19) 


Because  Es/N0  =  EB/N0  +  101og10£, 

Er/N0  =  Es/N0  +  101og107W -  101og10M.  (20) 

If  q  symbols  are  transmitted  contiguously,  then 
(18)-(20)  apply  when  TW  <  qM.  If  TW  >  qM,  then 
the  radiometer  parameters  are  TW  and  ER/N0  = 

^s/Nq  +  101og10#. 

When  the  symbol  error  probability  Ps  and  the 
number  of  bits  k  are  specified,  the  required  Es/N0 
for  CCSK  is  determined  by  solving  (4).  Equation 
(20)  is  then  solved  for  S  =  ER/N0  and  (16)  and  (17) 
are  used  to  determine  the  detection  probability  PD 
achieved  by  the  radiometer  for  a  given  PFA.  Results 
presented  later  take  into  account  the  effect  of  holding 
the  false-alarm  rate  (FAR)  a  constant. 

We  first  assume  that  the  radiometer  integrates  over 
Ts  seconds,  the  length  of  one  CCSK  symbol,  with 
the  integration  interval  matched  to  the  symbol.  This 
results  in  Es  =  ER  and  radiometer  time-bandwidth 
product  TW  =  M.  Table  II  shows  the  detection 
probability  PD  for  the  indicated  number  of  bits  k, 
when  Ps  =  0.001  and  PFA  =  0.0001. 

The  data  in  Fig.  10  exhibit  further  the  vulnerability 
of  CCSK  to  detection  by  a  radiometer  when  TW  =  M. 
Three  sets  of  curves  are  shown  for  three  choices  of 
PFA.  Within  each  set,  the  symbol  error  probability  Ps 
varies  from  0.00001  to  0.1.  From  these  sets  of  curves 
it  is  obvious  that  vulnerability  to  detection  decreases 
with  M,  but  increases  with  increasing  radiometer  PFA 
or  with  decreasing  symbol  error  probability  Ps.  Also 
included  is  the  spectral  efficiency  R/W  (in  bits  per 
second  per  hertz  of  bandwidth)  [8,  pp.  282-284], 
which  is  a  measure  of  performance  of  modulation 
methods.  For  this  case,  the  spectral  efficiency  is  given 

R/wPs&t:  (2D 

Note  that  the  vulnerability  to  detection  increases  with 
increasing  spectral  efficiency.  Also,  the  values  of  PFA 
used  in  Fig.  10  are  likely  to  be  larger  than  would  be 
used  in  practice.  However,  the  trend  is  obvious:  a 
further  decrease  in  PFA  will  lead  to  a  decrease  in  PD 
when  other  parameters  remain  the  same  as  used  in  the 
figure. 

We  now  assume  that  symbols  are  transmitted 
continuously  (and  are  contiguous).  An  assumption 
more  realistic  than  the  one  above  (TW  =  M )  is  that  the 
radiometer  integrates  over  multiple  symbols.  In  this 
case,  its  probability  of  detection  increases  compared 
with  the  single  symbol  case.  Conversely,  if  the 
radiometer  integration  time  is  less  than  the  length  of 
one  symbol,  the  probability  of  detection  is  decreased. 
This  is  illustrated  by  the  data  in  Fig.  11,  which  shows 
the  probability  of  detection  by  a  radiometer  as  its 
time-bandwidth  product  (i.e.,  integration  time)  varies. 
For  each  k  (k  =  10,  12,  14,  and  16)  results  are  shown 
for  a  range  of  values  of  Ps,  and  the  circles  represent 
the  results  given  in  Table  EL  Note  that  PD  increases 
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Fig.  10.  Spectral  efficiency  and  radiometer  detection  of  CCSK. 


Fig.  11.  PD  versus  TW  for  radiometer  detection  of  CCSK. 


TABLE  H 

Probability  of  Detection  by  Radiometer  with  Integration  Time  7^ 


Number  of 
bits,  k 

Required 
Eb/N0  (dB) 

Es/N0  (dB) 

10 

3.97 

13.97 

0.00170 

12 

3.56 

14.35 

0.00051 

14 

3.24 

14.70 

0.00024 

16 

2.97 

15.01 

0.00016 

considerably  with  increasing  TW  and  decreasing 
Ps  for  k  =  10  and  k  =  12.  However,  for  k  =  14  and 
k  =  16,  the  variation  of  PD  with  TW  and  Ps  is  small. 
Thus,  to  simultaneously  achieve  LPI  performance  (low 
PD)  for  small  values  of  Ps  requires  the  use  of  long 
sequences  with  lowered  spectral  efficiency  R/W. 


C.  Radiometer  Detection  of  TCCSK 

When  TCCSK  is  used,  the  equation  for  ER/N0  is 
given  by  (19)  with  M  replaced  by  MT: 

Er/N0  =  Eb/N0  +  moSl0TW  -  lOlo g]0MT  +  101og10 jfc. 

(22) 

We  first  assume  that  the  radiometer  integrates  over  the 
truncated  symbol  length;  i.e.,  TW  =  MT.  Results  are 
only  given  for  k  =  16  (M  =  65536),  and  illustrate  the 
effect  of  truncation  on  vulnerability  to  detection  by  a 
radiometer. 

Fig.  12  shows  PD  versus  the  truncated  symbol 
length  Mt  for  the  case  where  TW  =  MT.  Three  sets  of 
curves  are  shown  for  three  values  of  PFA.  Within  each 
set  the  symbol  error  probability  varies  from  0.00001 
to  0.1.  Also  shown  is  the  spectral  efficiency.  The 
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Fig.  12.  PD  and  spectral  efficiency  versus  truncated  symbol  length  MT. 


Fig.  13.  PD  versus  TW  for  radiometer  detection  of  TCCSK. 


results  in  Fig.  9  and  Fig.  1 1  show  that  the  detectability 
of  TCCSK  increases  slightly  compared  with  CCSK 
when  the  transmitted  symbol  lengths  are  the  same. 
Also,  TCCSK  slightly  increases  the  spectral  efficiency 
for  small  symbol  lengths.  This  occurs  because  the 
number  of  bits  per  symbol  is  a  constant  for  TCCSK 
(16  in  the  case  shown),  but  varies  for  CCSK.  Again, 
the  values  of  PFA  used  are  likely  to  be  larger  than  in 
practice,  but  the  trend  is  obvious. 

Fig.  13  shows  PD  versus  TW  for  M  =  216  and  for 
Mt  =  210,  212,  and  214  and  PFA  =  0.0001.  For  each 
choice  of  MT,  curves  are  shown  for  a  range  of  values 
of  Ps.  We  again  assume  that  symbols  are  transmitted 
continuously  (and  are  contiguous).  A  comparison  of 
Fig.  13  with  Fig.  11  shows  that,  for  the  same  symbol 


length  (i.e.,  Mr  =  2k),  the  detectability  of  TCCSK 
increases  considerably  compared  with  CCSK  for 
Mt  =  210,  but  only  slightly  for  longer  sequences. 

D.  FAR 

For  the  detection  results  presented  above,  the 
radiometer  PFA  is  held  constant.  It  is  usually  more 
realistic  to  hold  the  FAR  constant,  especially  when 
comparing  systems  with  different  integration  times. 
We  define  the  FAR  as  the  average  number  of  false 
alarms  per  second  (or  some  other  unit  of  time).  Its 
effect  on  detectability  is  discussed  in  detail  in  [7]  and 
is  summarized  here. 

We  assume  that  the  integrator  shown  in  Fig.  9 
operates  in  the  integrate-and-dump  mode.  That  is, 
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Fig.  14.  PD  versus  TW  for  CFAR  radiometer  detection  of  CCSK. 


the  integrator  output  is  sampled  every  T  seconds,  the 
integrator  is  reset  to  zero,  then  integrates  over  the 
next  T  seconds,  etc.  Thus,  a  decision  is  made  every 
T  seconds  and  the  false- alarm  rate  is  FAR  =  PFA /T.  If 
a  fixed  FAR  is  required  and  we  compare  two  systems 
with  integration  times  7j  and  T2  (7j  ^  Tz),  we  must 
use  a  different  PFA  for  each  system.  This  is  true 
even  if  the  bandwidths  of  the  two  systems  are  equal. 
Conversely,  if  the  integration  times  of  the  two  systems 
are  equal,  the  FAR  is  independent  of  their  bandwidths. 
However,  detection  performance  does  depend  on  their 
bandwidths  because  (16)  and  (17)  depend  on  the  TW 
product. 

Fig.  14  illustrates  the  effect  a  constant  FAR 
(CFAR)  has  on  radiometer  detection  of  CCSK  and 
should  be  compared  with  Fig.  1 1,  in  which  a  constant 
PFA  (CPFA)  is  assumed.  In  Fig.  14,  FAR  =  10~7 
and  is  expressed  in  units  of  Tc\  thus,  PFA  increases 
linearly  with  FAR.  The  results  for  both  cases  (CFAR 
and  CPFA)  are  identical  for  TW  =  1024;  however, 
for  TW  >  1024,  PD  is  larger  for  CFAR  and  for 
TW  <  1024,  PD  is  larger  for  CPFA.  These  situations 
occur  because  PD  increases  with  increasing  PFA 
(and  vice  versa)  when  TW  and  ER/N0  are  constant. 
Also,  because  PD  is  bounded  below  by  PFA  and  PFA 
increases  with  increasing  TW,  PD  approaches  1 .0 
asymptotically  as  TW  increases.  This  effect  is  evident 
in  the  results  shown  in  Fig.  14. 

Although  not  shown,  similar  results  are  obtained 
for  detection  of  TCCSK  by  a  CFAR  radiometer.  That 
is,  if  the  results  in  Fig.  13  for  CPFA  were  compared 
with  the  corresponding  results  for  CFAR,  the  same 
conclusions  made  just  above  would  apply. 

V.  CONCLUSIONS 

A  low  probability  of  intercept  (LPI) 
communication  technique  known  as  CCSK  has  been 


described  and  discussed.  This  technique  uses  cyclic 
(circular)  shifts  of  a  base  function  / (t)  to  modulate 
a  carrier.  The  function  f(t)  has  the  property  that  its 
cyclic  autocorrelation  has  a  distinct  peak  and  low 
sidelobes.  The  receiver  estimates  the  position  of  the 
peak  correlation  of  the  received  signal  plus  noise  with 
/ (t).  If  there  are  M  resolvable  positions  the  estimated 
position  represents  B  =  log 2M  bits. 

The  base  functions  considered  are  binary 
sequences  of  + 1  s  and  -Is,  which  results  in  biphase 
modulation  of  the  carrier.  Three  different  methods  for 
generating  these  sequences  were  evaluated  in  terms  of 
their  PMR.  The  first  is  an  MLS  of  length  M  =  2*  -  1, 
which  has  the  property  that  its  cyclic  autocorrelation 
has  a  peak  of  M  -  1  and  sidelobes  of  magnitude 
1.  Because  the  MLS  only  represents  k  —  1  bits,  an 
MMLS  was  considered.  The  MMLS  is  obtained  by 
appending  a  +1  or  -1  to  a  MLS,  thus  resulting  in 
a  sequence  of  length  2k  that  represents  k  bits;  Both 
the  MLS  and  MMLS  have  well-defined  structure,  and 
this  led  to  the  consideration  of  a  random  sequence  of 
length  2k.  Results  showed  no  significant  difference  in 
communication  performance  when  Comparing  the  use 
of  an  MMLS  with  a  random  sequence. 

Simulation  results  were  obtained  to  show  that  that 
the  performance  of  CCSK  in  terms  of  probability  of 
symbol  error  Ps  and  required  Es/N0  can  be  measured 
by  using  equations  that  apply  to  MOS,  Also,  it  was 
shown  that  the  receiver  signal  processing  for  CCSK 
is  simpler  and  easier  to  implement  than  for  MOS 
because  only  one  cyclic  correlation  is  computed  for 
CCSK. 

A  generic  radiometer  system  was  defined  and 
equations  for  evaluating  its  detection  performance 
were  given.  These  equations  were  applied  to  evaluate 
the  detection  of  CCSK  and  TCCSK.  Results  given 
for  various  combinations  of  parameters  show  that 
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to  achieve  LPI  performance,  transmitted  symbols  of 
length  greater  than  about  2 12  are  required.  The  use  of 
TCCSK  results  in  a  data  rate  higher  than  CCSK,  but 
with  a  penalty  of  higher  detectability.  The  tradeoffs 
between  detectability  and  data  rate  for  a  particular 
application  must  be  made  on  a  case-by-case  basis. 
However,  we  conclude  that  CCSK  and  TCCSK 
provide  LPI  capabilities  and  are  techniques  that  are 
simple  and  easy  to  implement. 
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